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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

COMPARISON OF THEORETICAL AND EXPERIMENTAL ZERO-LIFT 
DRAG-RISE CHARACTERISTICS OF WING -BODY -TAIL 
COMBINATIONS NEAR THE SPEED OF SOUND 
By George H. Holdaway 

SUMMARY 


The zero-lift drag rise at low supersonic speeds computed by 
linearized theory for several wing-body-tail combinations is compared 
with experimental data. The experimental data were obtained by the 
free-fall technique covering a Mach number range of M = 0.8 to M = l«l4. 
The procedures used in the theoretical calculations are outlined in 
detail. The results indicate that the computation method is capable of 
computing the drag rise at low supersonic speeds to an accuracy of about 
20 percent. 

Implications of the theory with respect to the selection of area 
distributions which will give a drag reduction over a range of Mach num- 
bers are examined. For cases where fuselage indentation is not feasible, 
calculations indicate that drag reductions may still be possible by 
increasing the fuselage volume in front of and behind the wing. 


INTRODUCTION 


Recent experimental results (refs. 1, 2, and 3) have demonstrated 
that the transonic zero-lift drag rise of wing-body combinations is 
primarily a function of the magnitude and rate of change of cross- 
sectional area along the longitudinal axis. This concept was utilized 
in the referenced tests to reduce the drag rise by indenting the body 
of a wing-body combination. These experimental results may be considered 
to be a qualitative verification of the linear theory as developed in 
references 4 through 8. 

The purpose of this report is to examine the quantitative relation- 
ship between the theory and ejqperimental data from free-fall tests of 
several wing-body combinations. In addition, the implications of the 


2 


MCA RM A53H17 


theory axe examined vith regard to modification of the area distribution 
of a configuration to keep the drag rise low over a selected range of 
transonic Mach numbers. 

This investigation utilized test data covering a Mach number range 
of M = 0.8 to M = l.l4 and a Reynolds number range of 2,500,000 to 
17,000,000, depending upon the wing mean aerodynamic chord of the con- 
figuration tested. 


SYMBOLS 


An 


C Dq 

^D 0 


dimensionless coefficients defining the magnitude of the 
harmonics of a Fourier sine series 

speed of sound, ft/sec 

drag at zero lift 

zero-lift drag coefficient, — 

^^w 

. _ _ , . , AD 0 

zero-lift drag-rise coefficient, — — 


ADo 

d 

l 

M 

N 

q. 

s 

s s 

$w 

V 

x 


zero-lift drag rise at transonic speeds, 11 

fuselage maximum, diameter, in. 

fuselage or "body length, in. 

Mach number, ~ 
a 

number of terms or harmonics used in the Fourier sine series 
dynamic pressure, i PV 2 , lb/sq ft 

projection of Sg on a plane perpendicular to x axis, sq in. 

cross-sectional areas formed by cutting the configurations 
with perpendicular or oblique planes X, sq in. 

wing area, sq ft 

free-stream velocity, ft/sec 

distance along the x axis measured from the midlength 
position, in. 

Cartesian coordinates as conventional body axes 
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6 angle between the z axis and the intersection of the cutting 

planes X with the yz plane 

p Mach angle, arc sin i 

M 

p mass density of air, slugs/cu ft 


9 


transformation of the length x 
radians 


to radians, arc cos — 7— , 

1/2 


X 


♦ 


S T 

s" 


a series of parallel cutting planes intersecting the x axis 
at the Mach angle P 

(At M = 1.0 these planes are perpendicular to the x axis; 
i.e., parallel to the yz plane.) 


angle "between the y axis and the intersection of the cutting 
planes X with the xy plane, arc tan (^/m 2 -1 cos 6 ^ 


first derivative of the projected cross-sectional area, 
second derivative of the projected cross-sectional area, 


dS 

dx 

d 2 S 

dx 2 


THEORY 

Development of the Theory 


Contributions to the development of the theory were made by the 
investigations reported in references 4 through 8. The calculated drag 
rise of slender bodies of revolution and the analogy between the drag- 
rise equations and the induced-drag equations for a wing are presented 
in reference 4. The development of equations for bodies of revolution 
with minimum drag rise is presented in reference 5 which also illustrates 
the application of a sine series to evaluate the theoretical drag rise. 

The theory was extended to arbitrary systems in reference 6, which showed 
how these arbitrary systems could be represented by equivalent bodies of 
revolution. This latter theory was applied to lifting surfaces in 
references 7 and. 8. 

The application (and limitations) of these methods to the theoreti- 
cal computation of the drag rise of wing-body combinations is presented 
in reference 9» The methods of reference 9 were used to make the theo- 
retical calculations for this report. The theory merely applies to the 
wave drag; any local shock or separation effects which might occur due 
to shape modification were not evaluated in the development of the theory. 
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Concepts Leading to the Drag Equation 


The derivation of the drag equation for a wing-body- tail combination 
is based on the theory that the configuration may be represented by a 
series of equivalent bodies of revolution. This theory is dependent on 
a simplified relationship between source strength and cross-sectional 
area. This relationship is used in planar wing and slender body of revo- 
lution problems. Specifically, as was pointed out in reference 9> the 
source strength is assumed to be proportional to the normal component of 
the stream velocity at the body surface. The theory also assumes that 
the configuration is of a conventional type with thin symmetrical airfoil 
surfaces and a high fineness ratio body. Further exceptions and limi- 
tations to the theory are given in reference 9- 

The development of equivalent bodies of revolution will be illus- 
trated by using the configuration shown in figure 1. The wing-body-tail 
combination is cut by a series of planes which always intersect the longi- 
tudinal axis at the Mach angle p. In other words, these planes are 
tangent to Mach cones. The plane identified in figure 1 as Xi repre- 
sents one plane of a series of parallel planes which cut the configuration 
along the entire longitudinal axis. Each plane of this series intercepts 
the yz plane in a line which forms the angle 0i = 0° with the z 
axis. Similarly, planes X 2 form the angles 0 2 with the z axis and 
the yz plane. For any one cutting plane of the series of planes 
X 2 = f (02, p) the oblique cross-sectional area is projected on a plane 
perpendicular to the x axis. This projected cross-sectional area is 
plotted as a function of x. The resulting plot may be considered as 
representing the longitudinal distribution of cross-sectional area S(x) 
of an equivalent body of revolution for the series of planes 
X 2 = f (0 2 , p) . For any one value of p this process is repeated for 
other values of 0 ranging from 0 = 0 to 0 = 2jt. However, if the con- 
figuration is symmetrical with respect to the xy and xz planes, then 
equivalent bodies for 0 from 0 to it only need be obtained. For 
bodies of revolution the area distribution is independent of 0. 

With these concepts and with the use of the simplified relationship 
between source strength and cross-sectional area, the equation for the 
zero-lift drag rise as a function of the rate of change of cross-sectional 
area can be derived from equation (46) of reference 7 and written as: 


£D 0 = — £1— J J' S M (xi) S"(x 2 )ln|xi - x 2 |dxidx 2 d0 

where xi and x 2 are two different locations along the x axis, 


S n (x) 


d 2 S( x) 
dx 2 


( 1 ) 


( 2 ) 
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This generalized equation can be simplified by solving the double 
integral of the functions of x through a Fourier sine series in the 
same manner as used in reference 5. 

J J' S"(xi) S"(x 2 )Zn|x a - x 2 |dx 1 dx 2 = ^ nA n 2 (3) 

n=i 

(The factor ikk appears because the body length, Z, has been defined in 
inches for this report.) 

Where 

N 

S-(x)=#Ui Y A n Sin nCp (4) 

dx 2 Z_. 

n=i 

and 

9 = arc cos ( s') 

1 / 2 ^ 

then 

r>Jt 

A n = ^ — / S’(x) sin n<pdcp (6) 

l Jq 

where the coefficients are a function of 0, since S’(x) is a function 
of 0. With this solution the simplified equation can now be written as 


AD 0 


1 pV 2 

144 2 



(7) 


The computing procedure followed in applying the foregoing equations 
and theory to the determination of the zero-lift drag rise is presented 
in the Appendix of this report. 


CONFIGURATIONS AND TESTS 


Plan-view sketches of the models tested, and also the axial distri- 
bution of cross-sectional area normal to the longitudinal axis, are shown 
in figure 2. The different configurations will be referred to as models 
A, B, C, and D as follows: 

Model A: aspect ratio 4 triangular wing with fuselage and tail 

Model B: aspect ratio 3 straight wing with fuselage and tail 
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Model C: aspect ratio 6, 45° sweptback wing with fuselage and tail 

Model D: fuselage and tail (consisting of two vertical and two 

horizontal surfaces) 


General geometric data for all the models are presented in table I, with 
greater detail given for model D in figure 3* The fuselage and tail were 
the same for all models. The fuselage ordinates from the 8- inch to the 
139. 4- inch station are given by the equation shown in figure 3* All- 
cross sections of the fuselage were circular, and the nose of the fuse- 
lage was faired from the ordinate at the 8-inch station down to a 1-1/2- 
inch-diameter boom 15 inches forward of that station. The fuselage was 
not indented. 


The free-fall technique employed is described in detail in refer- 
ence 10. The tests covered a Mach number range of 0.8 to 1.14 and had 
approximately the following Reynolds number variation: 


Model M = 0.8 M = 1.1 


A 

B 

C 

D 


7,500,000 

5.750.000 

2 . 600.000 

35,800,000 


16,600,000 

12,400,000 

5,600,000 

77,000,000 


The values of Reynolds number for models A, B, and C are based on the 
wing mean aerodynamic chords, and the values for model D are based on 
the fuselage length of 210.5 inches. 


The accuracy of the measurements of the drag rise divided by q 
and the Mach number is believed to be within ^2. = ±0.06 and M = ±0.01. 

q. 

Without an attempt to evaluate here the accuracy or limitations of line- 
arized theory, the theoretical computations as made are estimated to be 
accurate within 2 percent, except where linear theory might indicate a 
singularity in the drag curve. This accuracy of the theoretical data is 
based upon approximate methods for integrating the areas and determining 
the slopes of the area plots. Greater accuracies might be obtained by 
methods more exact than those outlined in the Appendix but probably are 
not justified. 


RESULTS 


The experimental zero-lift drag coefficients for the four models are 
presented in figure 4. The drag coefficients are based on the wing areas 
for models A, B, and C, and the fuselage maximum cross-sectional area for 
model D. The incremental drag rises above the subsonic values at M = 0.8 
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are presented in figure 5, together with the results of the theoretical 
computations. The rather large differences in wing area make comparison 
between configurations difficult; therefore, the data of figure 5 were 
multiplied by the respective areas used in computing the drag coef- 
ficients and the results presented in figure 6. The apparent late drag 
rise for model A is attributed to possible error in fairing the experi- 
mental data points which had scatter equal to = 0.06 from the mean. 


DISCUSSION 

Agreement Between Experiment and Theory 


Comparison of the drag-rise values presented in figure 6 indicates 
that, in general, the theoretical computations give a good prediction of 
the experimental values. Limitations in the theory resulted in only 
qualitative agreement being obtained between theory and experiment at 
Mach numbers below about 1.02. At Mach number 1.0 the experimental drag 
rise is less than the theoretical value, varying from 46 to 74 percent 
of theory. It may be noted that the greater the departure of the con- 
figuration from that of a body of revolution the greater the difference 
between theory and experiment became for M = 1.0. 

For Mach numbers greater than 1.0, the effect of theory limitations 
is diminished and the level of the drag rise is predicted quantitatively. 
In general, the values above M = 1.02 are in error by about 20 percent 
with a maximum deviation of experiment from theory of 26 percent. 


Application of Theory 


The prior section illustrated the adequacy of the theory in estimat- 
ing the zero-lift drag rise of wing-body- tail combinations. It is of 
interest next to examine the implications of the theory with regard to 
possible drag reductions through modifications of the area distribution. 

Consider a design problem involving a configuration similar to 
model B with an engine or other components within the fuselage which 
might make fuselage indentation impractical. In a case such as this the 
drag reduction would have to be attempted by adding volume before and 
behind the wing. That such an approach might be successful is illus- 
trated in figure 7(a). The solid curve is the computed drag rise for the 
unmodified configuration. The lowest, straight line represents the theo- 
retical drag rise for a Sears -Haack body with the same maximum cross- 
sectional area as the wing-body combination. The large difference 
between these two curves suggests the possibility of reducing the drag 
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rise of the original configuration by modifying its area distribution 
to approach, as a limit, the Sears -Haack distribution. 

The modification curves presented in figure 7(b) represent various 
fuselage revisions, obtained by adding volume to the fuselage in front 
of and behind the wing. Modification 1 is similar to the type of fuse- 
lage revision discussed in reference 1, in that it is based entirely on 
the distribution of cross-sectional area perpendicular to the x axis. 

Such a revision would be expected to reduce the drag rise at a Mach 
number of 1.0, and this anticipation is realized in figure 7(a). However, 
the desired equivalent body of revolution for M =1.0 is not the desired 
shape for other Mach numbers, and a rapid rise in drag with increasing 
Mach number therefore results as shown in figure 7(a). This drag rise is 
independent of possible separation effects which were not evaluated in 
the theory used. 

Clarification of the above point is afforded by consideration of 
figure 8(a). On this figure are shown the projected cross-sectional 
area distributions (S) of the original wing-fuselage-tail configuration 
for various angles *max- The largest angle *max is used for each 
Mach number computation and is the complement of the Mach angle. (.The 
equivalent tail area for the larger angles showed only slight differences 
from the tmax = 8.08°, so was omitted from this figure for clarification. 
The fuselage having a high fineness ratio had essentially the same area 
distribution for the entire Mach number range.) Shown on figure 8(a) is 
also a Sears-Haack body area distribution for a body with maximum cross- 
sectional area equal to that of the original wing, plus body. Modifi- 
cation 1, which optimizes the area distribution on the original wing- 
fuselage combination for a Mach number of 1.0, consists of additions to 
the body cross-sectional area sufficient to fill in the difference indi- 
cated by the shaded area. The fact that excess area is added for other 
Mach numbers is indicated by the penetration of the (s) curves for 
other values of 'I max into this shaded area. If it is desired to reduce 
the drag at Mach numbers somewhat above 1.0, much less volume should be 
added to the regions immediately in front of and behind the wing. 

In arriving at such a compromise for the wing body under consider 
ation, the following approach was followed to establish the lines for 
modification 2. The five area distributions (S) shown in figure 8(a) 
were arbitrarily replaced by a single curve representing the arithmetic 
average of the five curves. The Sears-Haack body area distribution was 
modified to represent a Sears-Haack body having a maximum cross-sectional 
area equal to the maximum of this averaged curve. The volume added to the 
fuselage was determined from the difference between this new Sears-Haack 
body curve and the average curve as shown in figure 8(b) . The ^max 
curves for M = 1.0 and M = l.l4 are also shown on figure 8(b) to illus- 
trate the reduction in the size of the new Sears-Haack body as well as 
the reduction in volume to be added immediately before and behind the 
wing* The computed drag rise for this second modification (fig* 7( a )) 
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shows an improvement in comparison to the original unmodified configur- 
ation over the entire design Mach number range of from 1.0 to l.l4. 

Modifications 3 and 4 were derived by the method described in 
reference 9 t° r the design of a fuselage modification optimized for a 
specific Mach number. Mach numbers of 1.05 and l.l4 were arbitrarily 
selected and the results of the calculations are presented in figures 7 (a) 
and 7(b). Modification 3 produced only a slight improvement over modi- 
fication 2, but it was maintained over the Mach number range shown. 
Apparently, as the minimum drag curve for a configuration is approached, 
slight streamwise modifications to the fuselage about the optimum shape 
can be made without greatly changing the drag. Further reduction of 
the drag rise at the higher Mach numbers at the expense of an increase 
at the lower Mach numbers could probably be effected by optimizing for 
the highest operating Mach number as was done in modification 4. Due 
to the limitations of the drag computations near Mach number 1.0, 
experimental investigations should be made to see if these small drag 
differences for modifications 3 and 4 really do exist. 


CONCLUSIONS 


The results of the computations and experimentation presented in 
this report indicate that for the models tested the following conclusions 
are justified: 

1. A computation method has been established which is capable of 
predicting the zero-lift drag rise at low supersonic speeds for a wide 
variety of wing-fuselage-tail combinations within an accuracy of about 
20 percent. 

2. In establishing the modifications to be made to a fuselage to 
reduce the drag rise, the effect of various revisions should be examined 
for a range of Mach numbers . 

3. For cases where fuselage indentation is not feasible, theory 
indicates that drag reductions may still be possible by increasing the 
fuselage volume in front of and behind the wing. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 17, 1953 
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APPENDIX 

PROCEDURE USED FOR COMPUTING ZERO-LIFT DRAG RISE 


The detailed application of the drag equation developed in the 
theoretical section of this report is outlined in the following para- 
graphs. The general procedures are presented with various simplifi- 
cations and suggested methods based on experience gained in making the 
computations for this report. The simplest case is presented first in 
order to better illustrate the basic steps and to shorten the presen- 
tation of the more involved procedures. 


Calculation at M = 1.0 for both Symmetrical and 
Unsymmetrical Configurations 


At M = 1.0 the cutting planes X are perpendicular to the x 
axis and the projected cross-sectional areas are independent of the 
angle 9. Therefore, the drag rise at M = 1.0 is computed in the same 
manner for symmetrical or unsymmetrical configurations. 

Step 1.- Determine the total cross-sectional area of the wing -body- 
tail combination and plot the area distribution as shown in figure 9(a)- 

Step 2.- Differentiate the area plot, figure 9(a) > with respect to 
x, and plot the results as a function of 9 as is indicated in 
figure 9(b) . The technique of reference 11 may be used to differentiate 
the area curve. 

Step 3.- Fit a Fourier sine series to the curve of figure 9(b) in 
order to determine the magnitude of the coefficients of the various har- 
monics. For example, from equation (6) of the theoretical section the 
coefficients A n are computed as follows: 


Ai 



9d 9 


Ao = — — 1^— sin 2 9d9, etc. 

« lJ Q dx 

The integrations may be accomplished by punch card machine computation, 
electronic-wave-analysis equipment or by ordinary long-hand computing 
methods. However, early investigations in connection with the compu- 
tat ions made for this report showed that a large number of terms (over 12) 
would be required to fit the irregular curves normally obtained, which 
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would rule out long-hand computing methods. Available electronic-wave- 
analysis equipment gave the relative magnitude of the coefficients but 
not their absolute values. The punch card computing methods were used 
for this report. Practical considerations of over-all accuracy did not 
justify calculating for more than 24 harmonics. On the other hand, com- 
puting for less than 12 harmonics reduced the resultant theoretical drags 
as much as 40 percent for some of the configurations investigated in 
this report. 

Step 4 .- Compute the drag rise from equation (7) which for M = 1.0 
reduces to the following coefficient form: 


Calculation for Mach Number Range, Configuration Symmetrical 
With Respect to xy and xz Planes, Wings 
and Tail Surfaces in xy Plane Only 

The first three steps in this calculation are essentially the same 
as for the M = 1.0 calculation except the computations should be made 
for at least four equivalent bodies of revolution rather than for just 
one. 

Step 1 .- Select the Mach number range to be covered by the compu- 
tations and determine the cross-sectional areas of the equivalent bodies 
of revolution for each value of 6 (p = constant) and plot as shown in 
figure 9(a) • By judicious selection of cutting planes for the compu- 
tation of the drag for the highest Mach number, the drag at lower Mach 
numbers can be computed from the same cross sections. This means that 
the values of Vy used for the lower Mach numbers are used for the 
intermediate values of ^ for the higher Mach numbers. For thin air- 
foil surfaces, which for purposes of this computation may be assumed to 
lie in the xy plane, the surfaces may be cut at angles ^ by planes 
perpendicular to the xy plane. The actual obliqueness or error in the 
wing cross sections disappears as the areas are projected onto the yz 
plane. It should be recalled that ^ is the intercept angle of a Q 
plane with the xy plane. The prior statements can be illustrated by 
an example showing the angles ^max which might be used in determining 
the cross sections for computations for M = 1.30. 



n=i 


In using equations (6) and (7) they may be simplified by deleting 2 /Z 
and (z/2) 2 , respectively, provided A n is considered as a dimension in 
inches . 
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max 


= arc tan //m 2 -1 


M 


^max> 


1.00 0 


1.01 8.08 

1.02 11.38 

1.05 17.75 

1.14 28.70 

1.30 39.71 


The equivalent bodies of revolution are determined for each of the above 
angles by cutting the configuration by a family of parallel planes. Then 
the area for each cross section is projected onto the yz plane by multi- 
plying the area by cos^rnax* This produces the equivalent bodies of 
revolution used in the final computations, where the drag computations 
for M = 1.05 will use only the equivalent bodies for the first four 
angles . 

For high fineness ratio fuselages, it is not necessary to cut the 
fuselage at various angles, because the projection of the cross-sectional 
areas onto the yz plane results in essentially a constant area distri- 
bution with Mach number as mentioned previously. For a fixed area distri- 
bution, the shape of the cross section of a high fineness ratio body is 
relatively unimportant if the body ends in a cylinder or a point (see 
ref. 12) . The zero-lift drag-rise coefficients for the Sears-Haack 
bodies can be computed from equation (l4), reference 5* which illustrates 
that the supersonic drag is constant: 


Step 2 .- Differentiate the area plots and plot as a function of 9 
as indicated in figure 9 (b) * 

Step 3 .- Fit a Fourier sine series to the curves of the area plots 
in the same manner as for the M = 1.0 computations. 

Step 4.- For the largest Mach number and angles compute the 
angles 0 and plot ZnA n as shown in figure 9(c). The same data can 
be used for the intermediate Mach numbers as follows: 
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M = l.l4 

M = 1.05 

M = 1.02 

A 

deg 

Gy 

radian 

ZnA n 2 

deg 

e, 

radian 

LnAjf 

J b 
deg 

e, 

radian 

ZnAftp 

0 

1.571 

90 

0 

1.571 

90 

0 

1.571 

90 

8.08 

1.285 

70 

8.08 

1.092 

70 

8.08 

.787 

70 

8.08 

1.857 s 

70 

8.08 

2.050 

70 

8.08 

2.355 

70 

II.38 

1.173 

60 

11.38 

.890 

60 

11.38 

0 

60 

11.38 

1.969 s 

60 

11.38 

2.252 

60 

11.38 

3-142 

60 

17.75 

.946 

50 

17.75 

0 

50 




17.75 

2.196 s 

50 

17.75 

3.142 

50- 




28.70 

0 

45 







28.70 

3 . l42 a 

45 








a Angles G are taken symmetrical about G = jr/2. 


Step 9 .- Integrate the area under curves similar to figure 9(c) and 
compute the drag-rise coefficient from the following equation derived 
from equation (7): 


n=24 


n=i 


d G 


Calculation for Mach Number Range, Configuration Symmetrical 
With Respect to xy and xz Planes, Wings 
and Tail Surfaces in Both Planes 


The essential differences from the prior method is that the verti- 
cal surfaces are rotated 90° into the xy plane, and if the same values 
of t used for the cutting of the horizontal surfaces are used to cut 
vertical surfaces then the cross sections for the vertical surfaces will 
correspond to different angles G than those used for the horizontal 
surfaces. Because the areas must be combined to give an equivalent body 
of revolution for one value of 0, cross plots of the areas should be 
made or the vertical surfaces should be cut at different angles than the 
horizontal surfaces. A satisfactory and simple procedure is to cut the 
vertical surfaces for one or two angles more than the horizontal surfaces 
to ensure a uniform variation of angles from 0 to 90°, and then proceed 
as follows: 

Step 1 . - Multiply the cross-sectional areas for the vertical 
surfaces by cost. 

Step 2 .- For each fuselage station plot the areas from step 1 as a 
function of 0y = - 90°- Subscripts V and H refer to the verti- 

cal and horizontal surfaces, respectively. 
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Step 3-- Read the areas from step 2 that correspond to the values 

of e~ 

Step 4.- Combine the cross-sectional areas (already multiplied by 
cos ty) for the vertical and horizontal surfaces to give one equivalent 
body for one angle Or = 9. 

Step 5 .- Compute the drag-rise coefficients as before for the case 
where all surfaces were in the xy plane. 


Calculation for Mach Number Range, Unsymmetrical Configuration 


This calculation is similar to the prior method except that the 
final integration must be performed from 6 = 0 to 9 = 2rt and the 
equivalent bodies of revolution determined from negative angles of t 
as well as from positive angles. These factors lengthen the problem 
but do not add to the complexity. 


Time to Perform Calculations 


Computations for this report, as was mentioned previously, were 
made by punch card machine computations . Time to compute the drag rise 
for a wing, similar to the models tested for this report, in combination 
with a high fineness ratio body from M = 1.0 to M = 1.2, was found to. 
require approximately 80 computer hours. This value is based on experi- 
ence gained in making the computations for this report and includes the 
time to lay out the cross sections and integrate the areas graphically. 
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TABLE I.- DIMENSIONS OF TEST MODELS 
(a) Wings 1 



Model A 

Model B 

Model C 

Area, sq ft 
Mean aerodynamic 

30.07 

21.68 

9.02 

chord, ft 

3.656 

2.813 

1.272 

Wing span, ft 

10.97 

8.20 

7.36 

Aspect ratio 

4 

3.1 

6 

Taper ratio 

0 

0.39 

0.5 

Airfoil section 

NACA 

0 to 0 . 5 c Ellipse 

NACA 64A010 


0005-63 

0.5 to 1.0c Biconvex 

perpendi cular 

Sweep of 0.25 

sStrearawise 

^-leading 

t/c =3 percent 

to 0.25 chord 

chord 

Twist, dihedral 

45° 

0° 

45 ° 

and incidence, deg 

0 

0 

0 


(b) Fuselage and tails for models A, B, C, and D 2 


Fuselage 


Fineness ratio 12.4 

Maximum diameter, ft 1.417 

Nose boom diameter, in 1*5 


Horizontal surfaces 


Area, sq ft 6.0 

Aspect ratio ^-5 

Taper ratio 0.20 

Airfoil section NACA 65 OO 6 streamwise 

Sweep of 0.25 chord ^5° 

Vertical surfaces 

Area, sq ft 3*3 

Aspect ratio 5*1 

Taper ratio 0.22 

Airfoil section NACA 65 OO 9 perpendicular to 0.25 chord 

Sweep of 0.25 chord ^5 


-‘-See figure 2. 
2 See figure 3- 



NACA RM A53H17 


17 




18459 


Figure 1.- Illustration of the cutting planes X and the angles 
9 and which are the intercepts of these planes with the 
yz and xy planes, respectively. 
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Figure 2.- Plan views of the models tested shown with their axial distributions of cross-sectional 

area normal to the longitudinal axis. 
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Figure 3* _ Geometry and dimensions of model D. 


vo 


MCA RM A53H17 


Zero-lift drag coefficient, Cp 


20 


NACA RM A53H17 



Mach number, M 

Figure 4 .- Experimental zero-lift drag coefficients, based upon wing 
area for models A, B, and C and upon maximum cross-sectional area 
for model D. 
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(a) Model A. 
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(d) Model D. 

Figure 5.- Comparison of experimental drag rises with values calculated 

utilizing linearized theory. 
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Figure 6 



(a) Model A. 
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(c) Model C. 
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(d) Model D. 

. - Comparison of theoretical and experimental increases in drag 
divided by dynamic pressure. 



Mach number, M 


(a) Calculated drag rise, Z£t, . 

v o 

Figure 7-- Comparison of the computed drag rise of model B with the computed values for the same 
configuration with fuselages of greater volume modified in various ways in an attempt to 
reduce the drag rise: 

(1) Enlarged body indented for equivalent wing area at M = 1.0 

(2) Enlarged body indented for the arithmetic average of equivalent wing areas (S) from 

M = 1.0 to M = 1.14 

(3) Enlarged body optimized for M = 1.05 

(4) Enlarged body optimized for M = l.l4 
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(t>) Fuselage revisions. 
Figure 7 • “ Concluded. 
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(a) Modification 1 for M = 1.0. 


Figure 8.- Modifications to the area distribution of model B to reduce drag by adding volume 

to the fuselage ahead of and behind wing. 
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(b) Modification 2 for M = 1.0 to 1.14. 
Figure 8.- Concluded. 
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0 = 1 arc cos radians 

373 ^ M 2 -i 

(c) 


Figure 9«~ The three steps in the theoretical solution of the drag rise. 
For an unsymmetrical configuration, plot C should he made from 
0 = 0 to 0 = 2rr. 
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